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Electronic screening renormalizes the linear bands of graphene and in the vicinity of the Dirac 
point, creates a diamond shaped structure in the quasiparticle spectral density. This is a result of 
electron-plasmon scattering processes which produce a finite momentum feature referred to as the 
'plasmaron ring'. In this work we explore the effects of uniaxial strain on these spectral features 
with the aim of understanding how strain modifies correlations. We derive and calculate the spectral 
density to the GoW-RPA level which allows us to identify the dispersive behaviour of the diamond 
geometry, and thus electron-plasmon scattering, for variation in electron-electron coupling strength 
and magnitude of applied strain. We find that the application of strain changes the geometry (in 
momentum) of the electron-plasmon scattering and that renormalizations beyond simple geometrical 
scalings further enhance this effect. These results suggest that the properties of the plasmaron ring 
can be tuned through the application of uniaxial strain, efl'ectively producing a larger fine structure 
constant without the need to change the sample substrate. 

PACS numbers: 79.60.-i, 72.10.-d, 73.21.-b, 73.22.Pr 



I. INTRODUCTION 

Graphene is a two-dimensional system comprised of a 
single layer of carbon in a honeycomb lattice. Near spe- 
cific points in reciprocal space, named K and K' valleys, 
the dispersion is conical, terminating at a Dirac point in 
each valley. Angle-resolved photoemission spectroscopy 
(ARPES) work by Bostwick et a\X has resulted in a sig- 
nificant interest in many-body interactions and their ef- 
fects near the Dirac point of graphene. In that work, 
and others^ it was observed that, for certain substrates, 
the single Dirac point splits into two pieces. In addition 
to this, the dispersion away from the K and K' points 
shows new structure due to electronic scattering with 
plasmons. This additional energy structure which occurs 
between the two Dirac-like points, illustrated schemati- 
cally in Fig.[TJa) at energy Ei, has been referred to in the 
literature as a plasmaron ring.— This structure has been 
shown to be well described by the GqW-RPA calculation 
of the screened electron-electron interaction. Further, 
there is some speculation that these additional energy 
structures, whose physical origin are a strong electron- 
plasmon scattering resonance, may be useful in the de- 
sign of plasmonic devices 

Recent infrared near-field optical techniques^"— have 
managed to observe plasmon dispersions and arc hoped 
to be complementary to existing ARPES work. Since 
these spectroscopic techniques are in general a momen- 
tum averaged probe, it is essential to have multiple tun- 
ing mechanisms in order to identify and separate the im- 
pacts of different many-body effects on spectral features 
as has been done in ARPES ii"— The most commonly ex- 
ploited tuning parameter is the chemical potential which 
can be controlled through gating as done, for example, 
in opticsi^iiS In the case of electron-electron interactions. 



the substrate also plays a substantial role, as it modifies 



the efltective fine structure constant, a 



through 



a change in the value of the effective dielectric constant 
of graphene at the substrate boundary, eo- Here g is 
the spin- valley degeneracy factor (g ~ gsgv = 4), e is 
the electron charge and u_f is the Fermi velocity. Since 
a is dependent upon the graphene substrate, tuning a 
is difficult in experiment, as one is somewhat limited in 
the values of a one can obtain from otherwise suitable 
substrates 



There is recent interest in strain on Dirac fermions 
in artificial honeycomb lattices produced in cold atom 
systems, which in some experimental set-ups naturally 
produce a strained Dirac dispersion.— In physical 
graphene, deformations as large as 20%i^"— are possi- 
ble above which an energy gap opens in the electronic 
spectrum. It may also be the case that strain can act 
as a tuning parameter for certain measurable quanti- 
ties and could be an important experimental tool. In 
this paper we derive and calculate the GqW-RPA re- 
sult for the screened electron-electron interaction (EEI) 
in doped graphene for arbitrary magnitude and direction 
of strain. We find that the application of strain modu- 
lates the radius of the plasmaron ring in fc-space and that 
the electron-electron interaction further enhances this ef- 
fect. These results suggest that the properties of the plas- 
maron ring can be tuned and therefore probed through 
the application of uniaxial strain, producing an effec- 
tively larger fine structure constant along the strained 
axis without the need to change the sample substrate. We 
begin with the basic formulas for a strained Dirac Hamil- 
tonian in Sec. [TTl In Sec. IIIII we discuss how the strain 
affects existing theory on the Gq-RPA level. Sec. llVlwill 
contain our main results, and Sec. |V] will conclude. 
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II. STRAIN 

The low energy Hamiltonian for a strained graphene 
sheet can be written in a hnear form for an effective mo- 
mentum, fe, as 



H = hvpc • k. 



(1) 



From this starting point, a system with strain appUed 
at an arbitrary angle 7 can be related to a symmetric 
system through the transformation^^ 



k = A(7)fc, 



(2) 



where ^(7) is the product of rotations and strain given 

by 



A(7) - i?(7)5(£)i?(-7)- 



(3) 



S{t) is a 2x2 straining matrix and i?(7) is a standard 
2-dimensional rotation through the angle 7. In the case 
of uniaxial strain, the strain matrix is a scaling in and 
ky given simply as 



5(e) = 







(4) 



where c\\ and c± are interrelated through some material 
specific constants: i"^ = %{^) — ^, related to the change 
to nearest-neighbour hopping, i, with a modification of 
the lattice spacing, a; Poisson's ratio v, the negative ratio 
of transverse to axial strain; and the strain modulus, e. 
We have taken these to have values of k = 1.1 and v = 
0.14 according to Ref. [l3 and arbitrary strain modulus. 
The relation to strain factors is given by 



C|| = 1 — 2Ke, 
cj^ = 1 + 2Kve. 



(5) 
(6) 



We illustrate schematically in Fig. [IJb) how a conical 
dispersion, with a circular cross section, when strained at 
an arbitrary angle 7 from the fc^j-direction produces an 
elliptical cross section. For any direction of strain, ^(7) 
is a 2x2 matrix, with elements Aij defined by the matrix 
multiplications of Eq. ([3|), which modifies the momenta 
k. One can then write an arbitrary vector k as 



fc 



All cos 9k + A12 sin 9k 
A21 cos 9k + A22 sin 9k 



Bx{9k, 
By {9k] 



(7) 



Here we have introduced a short hand for the x and y 
axis modifications due to strain, such that k — kB. For 
a given e and 7, the Aij are constant and one can then 
write a vector in terms of its strain contracted magnitude 
and direction. This allows for easy replacement of both 
the length and the angle of any unstrained vector with its 
contracted values. While in principle B contains strain 
modulus, strain angle and momentum angle dependence, 
for a given system we assume that e and 7 arc fixed, and 



(a) 



EEI 





FIG. 1. (Color online) (a) Inclusion of electron-electron in- 
teractions creates a diamond structure with energy width 5E 
and momentum radius, 5k, which describes the plasmaron 
ring, (b) Schematic illustrating a conical Dirac dispersion 
(blue) strained along an arbitrary angle 7 (red), with long 
and short axis factors cii and c 1 . 



so will refer to B as having only angular dependence. 
The length contraction is therefore simply 



B[9k)\^ jBl{9k) + Bl{9k) 



(8) 



where Bx and By are defined by Eq. ([7]) for a given an- 
gle 9k- The resulting angle of the strained vector k is 
therefore 



arctan 



By {9k ) 
Bx{9k) 



(9) 



Note that for improved clarity, we use the notation that 
barred variables have angles of /3 rather than 9. Simple 
manipulation would show that in the unstrained case, 
where C|| = = 1, you get the original magnitude and 
angle. Further, the momentum length contraction and 
angle shift are dependent upon the angle of the vector 
and of the applied strain. While these formulas and those 
that follow in Section IIIII are written for arbitrary mo- 
mentum and strain direction, it should be apparent that 
the cases of interest are for momentum either in the di- 
rection of the strain or transverse to the strain which we 
will refer to respectively as the long axis and the short 
axis. Depending on the application, one may wish to 
work purely in the principle axis coordinate system given 
by 7 = 0, which results in ^(7 = 0) = S{e). In this spe- 
cial reference frame, the strained momentum is 



(10) 



which is equivalent to the notation used elsewhereJ^ 





COS ( 






( c\\kx\ 




sin 


9k J 




\Cl_ky) 
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III. DYNAMICAL SCREENING AND SELF 
ENERGIES 

It has been shown in generaU^ that the density-density 
correlation function in the strained case, n(q,a;), is re- 
lated to its unstrained counterpart, n"(q,a;), through 

U{q,iu) = [detS{e)]-'U''{q,iu). (11) 

We see that the impact of strain is to both modulate the 
scattered momentum, q, and also scale the polarization 
by the inverse determinant of the strain matrix, S{e). 
Here, in the case of uniaxial strain, the determinant is 



simply C||Ci, representative of the modified slope of the 
density of states at the Fermi level. Further, since the 
polarization in the unstrained case only depends on the 
magnitude of q, the polarization in the strained case is 
evaluated at the scaled momentum given for a set of e 
and 7 in Eq. ([7]). 

The GoW-RPA self energy for the screened electron- 
electron interaction for a given momentum, fc, and fre- 
quency, Lo, in unstrained graphene has been previously 
derivedji^Ti^i The details of those derivations can be re- 
peated, but with the added complexity of strain which 
results in 



27r 



(12) 



and 



s'=±l 







^-^s.' (/?,,') 



2^ 



k+q 



^ fc+q fc+<J . 



(13) 



This calculation primarily samples magnitudes of scat- 
tered q over all directions dOq. The total self energy for 
a given band, s = ±1, is then the sum of these two con- 
tributions, E,(A;,a;) = Ef ^-^(fe, w) I]f''=(fc, w). With 
the exception of the additional strained labels these for- 
mulas [Eqns. ([T^ . p^ ] are identical to the unstrained 
self energiesJ^i^ Nevertheless we will outline precisely 
how these strain factors influence the evaluation of func- 
tions contained within Eqs. ()12m3p . To begin, the pref- 
actors are such that all quantities are scaled by the non- 
interacting chemical potential. Thus, the strength of the 
EEI rcnormalization increases or decreases with the value 
of the chemical potential. In the unstrained case, this 
scaling leads to momenta which are scaled by the Fermi 
momentum, kp- However, in the case of strain, since 
the slope of the dispersion is not uniform with angle, 
the Fermi momentum is instead given by the formula 
hkF\B{6k)\ ~ fi. Thus, in the strained case, hp depends 
on the angle of k and is shifted from the unstrained case 
by a factor of 1/15(6*^)1. Also note that the upper limit 
of the q integral requires a choice of cut-off to avoid an 
ultraviolet divergence. In principal one might be con- 
cerned with a variation of this cut-off for different direc- 
tions of k along the strained ellipse. The choice of cut-off 
has been shown to enter renormalization properties loga- 
rithmically, and for large cut-off these variations will not 
make a significant contribution to the result Despite 
this, we make the simplest correction for this error by 
modulating the cut-off by an angularly dependent factor 
of |-B(0fc)|. 



We can see the intrinsic dependence in the self en- 
ergy both on the direction of k and on the direction and 
strength of the applied strain. These factors enter into 
the strained dispersion, e| — , the inverse dielectric con- 

k+q 

stant, t~^{q,u})^ and the band overlap factors, Fss' ■ The 
inverse dielectric constant is given by 



'{q.^) = 



1 



q 



l-Vqliiq^u) q-all{q,ujy 
which through Eq. (|lip becomes 

q 



q - aU"{q,u})/{c\ic±)' 



(14) 



(15) 



where H" is the function in the unstrained case, but now 
evaluated at the strained momentum, q, derived repeat- 
edly elsewhere in explicit formi^i^i^ Vq = is the 
two-dimensional coulomb potential for a given effective 
dielectric constant of the medium, £q. 

The energies and wavefunctions are in the same form 
as in the unstrained case, but with strained variables 
for momentum and angle. The well known band overlap 
factor, Fss' , describing a transition from momentum k 
in the s band to momentum k' in the s' band acts to 
remove backscattering within a given band. One obtains 

^.s'fc) = ^[l + ss'cos(/3sfc-0] (16) 

where /3^^/ = — which can immediately be written 
down from Eq. 
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For clarity, we show an example of how (3/^1^/ coUapscs 
for the special cases of a fixed vector afong the long and 
short axes for the Of, ~ case. One requires two terms, 
the simplest is 



= arctan 



By 



arctan 



A21 
An 



(17) 



since A21 = for 7 = (long axis), and the same for 7 = 
^ (short axis). One could instead provide this example 
in the coordinates of the strain, as in Eq. (jlOp. which 
fixes 7 = 0. In this case, one finds 



Ik- 



arctan 



(arctan (^fi^ 



0, long axis, 9k = 
^, short axis, 9k = ^ 



(18) 



a different result from Eq. \T7l but in both cases one ob- 
tains in Eqn. ([TT]) and ([T5)) that = 9^, or that the 
strain transformation does not modify the vector angle 
for a momentum either in the long or short axial direc- 
tions. 

The k' term is more complicated, as it depends on 
the angle to the vector k' . Continuing our example for 
9k = 0, one would obtain 



9k' = arctan 



q sin 9q 
k + q cos 9q 



arctan (i?) (19) 



where we again have introduced a shorthand notation. 
We can then write the desired angle 



f3-k' 



arctan 



A21+A22R 
An +A12R 



(20) 



Finally, in the case of no strain, A21 = A12 = and 
An ~ A22 = 1 and wc find that 9k'- In this 

limit the overlap function Fgs'iPkk') ~^ Pss'{9kk') which 
is the standard result for unstrained graphenej^i^ We 
can see from the above examples that one can align the 
coordinate system along either the direction of strain, 7, 
or along the momentum in question, k. All that matters 
for the calculation is the direction of k relative to the 
strain. In practice, the two most interesting cases are for 
a momentum along the direction of strain (long axis) or 
perpendicular to it (short axis). 

The energy dispersion for a band, s, is of the general 
form 



s\k + q| - 1 



(21) 



which is in units of the bare chemical potential, /x. This 
can be written out in terms of fc, q and the strain angle 
by using Eq. ([7]) . We can then write the strained Green's 
function for each band 



G'{k,uj) 



cj - e|- - Es(A;,a;) 



(22) 



The total spectral function is A(k,uj) = ^ A^{k^Lij) 

s=±l 

which follows from Eq. (|22p and is 



A{k,uj) 



s=± 



1 



— lml]s(fc, uj) 



[lo - RcS,(fc,a;) - e|]2 + [Iml],(fc, uj)f ' 

(23) 

Without self energy, the spectral density in Eq. re- 
duces to a sum of two Dirac delta functions, and the 
strained dispersion curves are a simple geometrical scal- 
ing of the momentum axis applied to the bare unstrained 
dispersions. When one has interactions, or finite self en- 
ergies, as in Eq. (^5]) . one no longer obtains the strained 
spectral density by a simple geometric scaling of momen- 
tum of the unstrained spectral function. Of key impor- 
tance in the discussion that follows is that Eq. gives 
results which differ from simply taking the spectral den- 
sity in the unstrained case, which includes EEIs, and then 
scaling the momentum by an appropriate strain factor of 
\B\. Thus, we will use this simple geometrical scaling 
as a comparison to accentuate the impact of strain on 
correlation effects. 



IV. RESULTS 

Shown in Fig. [5] is an array of cuts along each of the 
long and short axis of the spectral function of a strained 
graphene sheet for the 10% and 20% strained cases for 
an EEl coupling strength of a = 2 (a value relevant for 
hydrogen passivated silicon carbide, H-SiC)ii We begin 
our discussion by considering the slice in energy through 
the A; = momentum point. In a non- interacting sys- 
tem, one would expect that the spectral function would 
exhibit no modification at the Dirac point from the ap- 
plication of strain. However, with interactions, the ap- 
plication of strain modifies the local electronic structure 
which produces additional renormalizations not due to 
the simple geometry of the strain. We show a specific 
interacting case in Figl^I^a). While there is no directional 
dependence at fc = which would distinguish the long 
and the short axial directions, there is nevertheless a no- 
ticeable strain dependence. Interestingly there is little 
to no modification of the primary Dirac point identified 
by cj = i?o. At larger negative energies the peak located 
at E2 shifts significantly, resulting in a 5E which grows 
under the application of strain. 

The full calculation along the long and short axis are 
shown in Fig.[21Jb-e). One observes that the renormaliza- 
tions do not influence the location of the Fermi level at 
a; = (shown by the white dashed line) beyond the sim- 
ple geometric modification along the momentum axis as 
discussed in Sec. lIIII At higher energies there are distinct 
features which go beyond a simple scaling of fc. To help 
see this, we take a characteristic slice of the colour plots 
through the energy E\ which we define as the midpoint 
between Eq and E2 ■ We then have a robust way of iden- 
tifying the radius of the plasmaron ring, Sk, and how it 
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FIG. 2. (Color online) {a,)A{k — 0, uj) for the long and short 
axis for increasing strain. SE is robustly defined as the energy 
separation between the two dominant peaks. The shift in E2 
would not occur in a non-interacting system, (b— >e) Color 
plots representing for cuts along the long and short 

axis for e = 0.1 and 0.2. (f) Complementary to the color 
plots: A cut through the diamond structure at u — Ei = 
{Eq + £2)72. The shift in 5k along the long axis is labelled. 



disperses with strain, e, and the EEI coupling strength 
a. These slices are shown in Fig. [5Jf). The application of 
strain results in an increase in Sk along the long axis (left 
frame) and little change for the short axis (right frame). 
The later effect in the short axis is simply a manifesta- 
tion of the relative values of C|| and c±. Along the long 




_i I I I I I I I I I li I I I I I I I I I lJ q 

0.5 1 1.5 2 2.50 0.5 1 1.5 2 2.5 

a a 



FIG. 3. (Color online) Left frame: Peaks in A{k = 0, cj) as a 
function of EEI coupling strength for strains of e = 0.0, 0.1, 
and 0.2. Peaks are identified as Eq, E\ and £2- Right frame: 
The energy width of the EEI diamond structure, defined as 
SE = E2- Eo. 



axis, one immediately sees features which do not allow 
the various strained cases to scale onto each other as 
would be the case without the EEI renormalizations. As 
e increases from 0.1 to 0.2, 6k increases and there is ad- 
ditional broadening of the peaks. Also, there is a change 
in the value of the spectral function along fc = which 
is important to note, since any modification at A: = is 
significant as it relates to changes in correlation effects. 

Continuing the discussion of the k = variation, we 
focus on the dispersive nature of the two independent en- 
ergy features, Eq and E2, which are identified as peaks 
in A{k = 0,w) as shown for an example in Fig. [51Ja). In 
Fig. [3] we show our results for these extracted peaks as a 
function of the EEI coupling strength, a, as well as for 
variation in strain. There is no difference between the 
long and short axis properties for this fc = cut. As 
noted, we find that strain shifts E2 to larger negative en- 
ergies while Eq remains essentially unchanged. This is 
interesting in that we can identify the upper cone which 
extends from fc = at a; = £"0 to /c = /cf at a; = in 
which, even with the strain-corrected self energy, we see 
only the simple modification due to strain. This is not so 
for other regions of k and uj, near E2- Physically this rep- 
resents how different regions^^^ of the inverse dielectric 
function are sampled in Eqs. (|12j|13p and suggests that E2 
is a feature dominated by coupling to the modified plas- 
maron band of Eq. (fTTI) while Eq (and the band extend- 
ing towards kp) is primarily due to screened electron- 
electron coupling in the intra or interband regions^^ and 
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FIG. 4. (Color online) Extracted values of Sk along the long 
axis as a function of EEI coupling strength, a, for unstrained 
e = 0.0, and strained graphene, e = 0.1 and 0.2. 5k is identi- 
fied as in Fig. [2jf). Dashed lines are 5k measured in the sim- 
ple geometric scaling case for scaling factors relevant to strain 
values of e = 0.1 and 0.2 (red and blue dashed respectively). 
The shaded regions accentuate the additional contributions 
to 5k from the correlation effects due to strain. 

therefore is modified in a simple way by the application 
of strain. The dispersive nature of E2 results in a depen- 
dence of SE on the strain amplitude, shown in the right 
hand frame of Fig. [3l which for small e is expected to be 
a linear change from the unstrained casei^ 

Next we address the momentum scaling of the plas- 
maron ring in more detail. We plot in FigU] how Sk is 
modified along the long axis with applied strain mod- 
ulus, e, and with variation in EEI coupling strength, 
a. We restrict a to a range of values where plasmaron 
bands have been observed expcrimcntallj*^ and empha- 
size that larger values of a do not produce as clear fea- 
tures of the plasmaron band due to increased scatter- 
ing as explored in previous calculations on unstrained 
graphenci ^^i^^ Shown in black for the unstrained case, 
this curve agrees precisely with previous calculations by 
Bostwick et aii which were found to agree well with ex- 
perimental results. Of primary interest in this work is 
the impact of strain, shown for e = 0.1 and 0.2 (red and 



blue solid lines respectively) across the same range of a. 
We compare these to the same feature in the case of a 
simple geometric scaling of momentum which introduces 
a 1/1 J5| = l/c|| scaling in Sk shown as the dashed curves. 
We highlight the difference between the simple scaling 
result and the full GqW-RPA calculation of the peak to 
peak distance Sk. The impact of the strained electron- 
electron interaction is to push the value of Sk to an even 
larger value than with momentum scaling alone. Finally, 
we mark an example (open circles) where the application 
of 10% strain to graphene at an EEI coupling strength of 
a = 1 produces a plasmaron ring with an effective Sk/ ^ 
equivalent to an unstrained sheet of graphene with an a 
value nearly 70% larger. 

V. CONCLUSIONS 

We have studied the dependence of the electron- 
plasmon scattering structures of the single particle spec- 
tral density, for variation in EEI coupling strength with 
a uniaxially strained Dirac dispersion at finite chemical 
potential. We have derived and calculated the GqW- 
RPA self energy and full spectral density under strain 
and compared to the case of a simple geometric scaling 
of momenta as would be valid in a non-interacting sys- 
tem. We find that along the long axis, the strained EEI 
acts to modify the Dirac point separation energy, SE, and 
also to create an enhancement of the geometrical strain 
modification of the primary feature of electron-plasmon 
coupling, the radius of the 'plasmaron ring', Sk. In addi- 
tion, we find that these extra strain-induced renormaliza- 
tions have very little effect on quasiparticle peaks close 
to the Fermi level. By probing the regions in k and lu 
where strain has an effect on the correlations in the spec- 
tral density, it should be possible to isolate the changes 
in many-body effects from those due solely to geometry 
and sample, such as substrate effects which maybe be 
of concern in ARPES experiments. Finally, we suggest 
that, in so far as the plasmaron ring is concerned, the 
application of strain could be used to tune the effective 
fine structure constant along the axis of applied strain in 
lieu of changing the sample substrate. This concept will 
be very useful for any momentum sensitive probe where 
the signal in the direction of strain can be isolated, such 
as ARPES. 
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